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Abstract 

Using a very high precision spectral calculation applied to the incompressible and 
inviscid flow with initial condition ^)q{x\, X2) = cosxi + cos 2x2, we find that the 
width S(t) of its analyticity strip follows a ln(l/i) law at short times over eight 
decades. The asymptotic equation governing the structure of spatial complex-space 
singularities at short times (Frisch, Matsumoto and Bee 2003, J. Stat. Phys. 113, 
761-781) is solved by a high-precision expansion method. Strong numerical evidence 
is obtained that singularities have infinite vorticity and lie on a complex manifold 
which is constructed explicitly as an envelope of analyticity disks. 
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1 Introduction 



In early September 2001 one of the authors (UF) attended the Zakopane meeting on Tubes, 
Sheets and Singularities (Bajer and Moffat, 2003) which was also attended by Richard Pelz. 
There were many discussions about the issue of finite-time blowup for 3D incompressible 
Euler flow. Richard, who had studied a flow introduced by Kida (1985), had evidence in 
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favor of blowup but one could not rule out that the highly special structure of this flow 
would lead to quasi-singular intermediate asymptotics. The three authors of this paper 
then decided to embark in a long-term project aimed at getting strong evidence for or 
against blowup for a wide class of flows encompassing the Taylor-Green flow (Brachet 
et al, 1983) and the Kida-Pelz vortex (Kida, 1985; Pelz, 1997; Pelz and Gulak, 1997; 
Pelz, 2001), namely space-periodic flow with or without symmetry having initially only a 
few Fourier harmonics. Such initial flow is not only analytic but also entire: there is no 
singularity at finite distance in the whole complex spatial domain. 

As has been known since the mathematical work of Bardos, Benachour and Zerner (1976), 
any real finite-time singularity is preceded by complex-space singularities approaching the 
real domain and which can be detected and traced using Fourier methods (Sulem, Sulem 
and Frisch, 1983; Frisch, Matsumoto and Bee, 2003). This method is traditionally carried 
out by spectral simulations which run out of steam when the distance 5(t) from the real 
domain to the nearest complex-space singularity is about two meshes. We pointed out in 
a recent paper (Frisch, Matsumoto and Bee, 2003, henceforth referred to as FMB), which 
also reviews the issue of blowup, that it may be possible to extend the method of tracing of 
complex singularities by performing a holomorphic transformation mapping singularities 
away from the real domain and, perhaps doing this recursively. This is the basic idea of 
the spectral adaptive method which aims at combining the extreme accuracy of spectral 
methods with the local mesh refinement permitted by adaptive methods. 

In one dimension the complex-space singularities of PDE's are isolated points, at least in the 
simplest cases, as for the Burgers equation. In higher dimension they are extended objects, 
such as complex manifolds. Understanding the nature and the geometry of such singularities 
is a prerequisite for mapping them away. Many aspects can already be investigated in the 
two-dimensional case for which we know not only that blowup is ruled out, but we also know 
that the flow stays a lot more regular than predicted by rigorous lower bounds (basically 
5(t) seems to decrease exponentially whereas the bound is a double exponential in time). 
In FMB we gave some evidence that in 2D the complex singularities are on a smooth 
manifold, but the nature of the singularities was not very clear and in particular the issue 
of finiteness vs. blowup of the complex vorticity was moot. 

In FMB we also pointed out that the issue of singularities can be simplified if we limit 
ourselves to short times. Let us briefly recall the setting. We start with the 2-D Euler 
equation written in stream function formulation 



where J(f,g) = d 1 fd 2 g — digd 2 f- As in FMB, we focus on the two-mode initial condition 



d t V 2 4> = J(^,VV), 



(1) 




(2) 



2 



one of the simplest initial condition having nontrivial Eulerian dynamics. The solution 
ip(z, t), obtained by analytic continuation to complex locations z = x + iy, is expected to 
have singularities at large imaginary values when t is small. If one focuses on the quadrant 
yi — > +00 and y 2 — > +00, an asymptotic argument given in FMB suggests looking at 
solutions satisfying the similarity ansatz 

*l>(z,t) = (l/t)F(z), (3) 
z = (zi, z 2 ) = (z 1 + i\nt, z 2 + (i/2)\nt). (4) 

Substitution in (1) gives the similarity equation 

V 2 (-l + id 1 + {i/2)d 2 )F = J(F, V 2 F), (5) 



where the overscript tilde means that the partial derivatives are taken with respect to the 
new variables. The initial condition (2) becomes a boundary condition 

F(£)~i(e- lil +e- 2ii2 ), y x -> -00, y 2 -> -00. (6) 



Note that (5) has no time variable. Its solution can be shown to be analytic for sufficiently 
negative y~\ and y 2 . If it has singularities in the real or complex domain then, by (4), the 
solution of the original Euler equation should have short-time singularities at a distance 
5(t) oc ln(l/t). 1 

The outline of the present paper is as follows. In Section 2 we check the logarithmic law for 
5(t) and thus the validity of the similarity ansatz. In Section 3 we develop a new technique 
for solving the similarity equation (5) in both the real and complex domains. In Section 4 we 
present numerical results on the nature and the geometry of the singularities. In Section 5 
we show how to actually construct the singular manifold from the Fourier transform of the 
solution. In Section 6 we make some concluding remarks. 



2 Numerical validation of the short-time behavior 



The standard way of measuring the distance 5 of the nearest complex-space singularity 
(also called the width of the analyticity strip) is to use the method of tracing (Sulem, 
Sulem and Frisch, 1983). Indeed, the spatial Fourier transform vj, of a periodic function 

1 This law, for flow having initially a finite number of Fourier modes, was first derived for the 
one-dimensional Burgers equation and conjectured to apply also to the multi-dimensional incom- 
pressible Euler equation (Frisch, 1984). 
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v{x) has its modulus decreasing at large wavenumbers k = \k\ as e~ Sk in both one and 
several space dimensions. More precisely, for each direction k = k/k, the modulus of the 
Fourier transform decreases as e _<5 fe fe ; the width of the analyticity strip is then the minimum 
over all directions k of 8^ which - by a steepest descent argument - also controls the high-/c 
decrease of the angle average of the modulus of the Fourier transform. 2 

Accurate measurement of the width of the analyticity strip 8 gets difficult when it becomes 
smaller than a few meshes, so that there is not enough resolution to see long exponential 
tails (Brachet et al, 1983). A different difficulty appears when 8 is very large and the 
Fourier transform decreases so rapidly that it gets lost in roundoff noise at rather small 
A;'s. To check on the logarithmic law of variation of 8(t) at short times for the full time- 
dependent Euler equation (1), we have to face the latter difficulty. To overcome it, we 
employed a 90-digit multiprecision spectral calculation 3 with the number of grid points 
ranging from 64 2 to 128 2 The temporal scheme is fourth order Runge-Kutta. Fig. 1 shows 
the k dependence at various short times of the angular averages 4 of the modulus of the 
Fourier transform of the velocity v = (d2ip, —diifj), where di and 82 are the derivatives with 
respect to x\ and Xi- Note that there is a strong odd-even wavenumber oscillation. This has 
to do with the interference of two complex singularities separated by 7r in the xi direction 
(a consequence of a symmetry of the initial condition). A very clean exponential decrease 




Fig. 1. Angular averages of the modulus of the velocity Fourier modes at short times calculated 
with 90-digit precision. Strong oscillations between odd and even /c's are due to the symmetry of 
the initial condition (2). The values of the time steps for integrating the Euler equation with 64 2 
and 128 2 grid points are respectively 2 x 10 -11 and 5 x 1CT 6 . 

is observed for even wavenumbers. This allows the measurement of 8(t) for hundreds of 
values of t covering a very wide range, as shown in Fig. 2 in log-lin coordinates. It is seen 

2 For the way 5^ is related to the singular manifold in the two-dimensional case, see Section 5. 

3 All multiprecision calculations in this paper were done using the package MPFUN90 (Bailey, 
1995). 

4 More precisely, we use "shell-summed averages" , defined at the beginning of Section 4. 
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Fig. 2. Width of the analyticity strip 5(t) at short times 2 x 1CT 10 < t < 7.3 x 10 _1 measured 
from the angular average of the modulus of the velocity Fourier modes restricted to even /c's. 

that, up to approximately t = 1CT 2 , the logarithmic law is satisfied over eight decades, a 
range which could not be achieved without the use of very high-precision computations. 



3 Solving the similarity equation 

We are interested in space-periodic solutions to (5) satisfying the boundary condition (6). 
This problem can be solved exactly using Fourier series 

oo oo 

F ( z ) = E E h lM e t{klzl+k2Z2) - (7) 

fci =— oo k2 =— oo 

Note that we have dropped all tildes on the space variables since from now on we will work 
exclusively with the similarity variables. Obviously, the boundary condition (6) allows the 
presence only of Fourier harmonics with wavevectors in the quadrant (k± < 0, k 2 < 0). 
After (5) is rewritten in terms of the Fourier coefficients of F, the multiplications appearing 
in the Jacobian go over into simple convolutions with only finitely many terms, because all 
wavevectors involved must have non-positive components. It follows that, for integer n > 1, 
all Fourier coefficients for wavevectors which are on the line 2k\ + k 2 = — 2(n + 1) (see 
Fig. 3) are expressible in terms of the finitely many Fourier coefficients with wavevectors 
lying above this line. More precisely, for any point P on this line, the region of dependence 
is the grey-shaded rectangle. Specifically, we define the selective Fourier sum over this line 

F n (z) = EF„ CT ,_ 2(n+1 _ (j)e -^ +2 ^ 1 - CT ) 22 ), (8) 

CT=0 
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-2(n+l) 

Fig. 3. Region in the wavenumber space relevant in calculating the Fourier coefficients Fk lt k 2 on 
the line 2k\ + k2 = — 2(n + 1). The largest disk centered at the origin fitting into this region is 
also shown. 

so that 

oo 

F(z) = 5> B (z). (9) 

n=0 

From (5) we then obtain the recursion relations 
1 n 

V 2 F n+1 = — ^ J(F r ,V 2 F n _ r ), (10) 
n + 1 r = 



with 



-IZl 



+ e 



-2iz 2 \ 



(11) 



Equation (10) can then be rewritten as a set of algebraic equations for the Fourier coeffi- 
cients: 



oo 

Fki,k2~ ^fc(p) F pij p2 -ffei— p!, fe 2 — P2 J (12) 



2 


k — p 


\ 2 (pAk) 


2fci + A; 2 + 2 


k 


2 



(13) 



where k A k' = kik' 2 — fe&i- One could solve the Poisson equations (10) which recursively 
define the F n 's using FFT methods as in Brachet et al. (1983). Alternatively - and this is 
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the method used here - equations (12)-(13) can be used to calculate recursively the exact 
expressions of all the Fourier coefficients. The first two F n 's are 



F 1 (z) = ~e< n+2n \ (14) 

F 2 (z) = - A e -^l+4*2) + ^_ e -i(2 21 +2 22 )_ ^ 

We have numerically determined the F n 's for n up to n max = 1500, using quadruple- 
precision (35-digit) accuracy. As we shall see in the next section, lower accuracy would give 
spurious results. From (12) and (13), the F kl ^ 2 are obviously real. Furthermore, we find 
numerically that their signs alternate: specifically (— l) hl F kl k2 > 0, for all (ki,k 2 ) except 
(-1,0). 



4 Numerical results on singularities 



We shall work here mostly with shell-summed (Fourier) amplitudes. For a given periodic 
function /, we define its shell-summed amplitude as 

K<\k\<K+l 



where the fk are the Fourier coefficients of /. The shell-summed amplitudes for the solution 
F to the similarity equation (5) can be calculated, in principle exactly, for all wavenumbers 
K such that K + 1 < k max where k mSLX = (2/y / 5)(n max + 1), this being the radius of the 
largest disk centered at the origin and contained in the triangular region of Fig. 3. 

Fig. 4 shows, as a function of K, the shell-summed amplitude of F for n max = 1500. As seen 
in Fig. 4(a), the results obtained with 15 and 35 digit precisions differ markedly beyond 
wavenumber 800. Even the 35-digit calculation becomes unreliable beyond wavenumber 
1300. In FMB, when we discussed results about singularities without resorting to short- 
time asymptotics, we reported various difficulties: the need to perform a kind of Krasny 
filtering (Krasny, 1986; Majda and Bertozzi, 2001; Caflisch, Hou and Lowengrub, 1999) 
and our failure to improve the range of scaling by going to resolutions higher than 512 2 . It 
is now clear that the key to clean results is to use high precision. 

We fit the shell-summed amplitude of F by a function of the form CK~ a exp(— SK) and 
find the best fit to be 0.5-fT -2 ' 16 exp(— 0.0065-fT), as shown in Fig. 4(b). The fit is done by 
least squares in the lin-log representation over the interval < K < 1100. 
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Fig. 4. (a) Shell summed amplitudes A K for the stream function F calculated with 15 and 
35-digit precisions in lin-log coordinates. The strange behavior of the 15-digit result in the high 
wavenumber region k > 800 is due to insufficient precision, (b) Shell summed amplitude A K 
calculated with 35-digit precision and almost overlapping least square fit to CK~ a exp(— 5K) 
in lin-log coordinates. Numerical instability is observed for K > 1300. Inset: with 35-digit 
precision and its least square fit in log-log coordinates. 

Hence, the solution to the similarity equation has complex-space singularities, the closest 
one being within 5 ~ 0.0065. This relatively small value, for an equation in which all the 
coefficients are order unity, is accidental and can be changed by slightly modifying the 
coefficients in front of the two harmonics in the initial condition (2) and in (11), as will be 
shown later in this section. 

The angular dependence of the mode amplitude at high wavenumbers can also be obtained. 
Fig. 5 shows that |Ffe| decreases exponentially with k = \k\ for a given direction k = 
(cos 9, sin#). 5 We see that the logarithmic decrement, denoted by 5(9), varies strongly 
with 9. Its variation is shown in Fig. 6 for n < 9 < 3n/2, that is over the third angular 
quadrant where k\ and k 2 are negative. 6 Obviously, the width of the analyticity strip is 




1400 



5 = min 6(9), 



;i7) 



the minimum being achieved in the most singular direction 9 = 9*, such that tan^ « 1.62. 
The function 6(9) can be related to the singular manifold (see Section 5). As in one di- 
mension, the prefactor of the exponential in the shell-averaged amplitude A K contains 
information about the nature of the singularities. In the inset of Fig. 4 we see about one 
decade of power-law scaling before the exponential falloff. This range can be increased by 

5 Actually, we sum over all Fourier modes for which k is within 7r/100 of the direction 0. 

6 At the edges of this quadrant 5(0) becomes infinite (W. Pauls, private communication), but 
this is hidden by the slight angular averaging. 
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Fig. 5. Amplitudes of the Fourier modes \Fk\ along five different directions from 1.057T to 1.327T 
(from bottom to top) 

10 
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0.001 
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e/n 

Fig. 6. Logarithmic decrement 5(6) as a function of the polar angle 6 in the (k±, fo) plane. The 
minimum of 5(6) is around 0.0065. The most singular direction is 6+ = 1.324-7T. 

moving closer to the singularity, more precisely to the point z+ on the singular manifold clos- 
est to the real domain which has imaginary part Imz* = (yi+, ?/2*) with yi± = 5 cos 0* and 
1/2* = <5sin#*. Such an imaginary shift produces a function F^ 1 ) whose Fourier coefficients 
are obtained by multiplying the Fourier coefficients of F by exp (hk\ cos 9* + hkz sin #*) , 
where < h < 5. 7 Fig. 7 shows shell-summed amplitudes of F^ for four values of h. For 
h = 5 more than two decades of power-law scaling is seen with an exponent somewhere 
between —2.0 and —2.2. This clean scaling evidence has an important consequence for the 
vorticity on the singular manifold (see below). 



7 Note that, in the short-time asymptotics, an imaginary shift (hi, /12) amounts to changing the 
initial condition (2) into ipo(z) = e hl cos(zi) + e 2fl2 cos(2z2) within terms irrelevant for yi — > +00 
and j/2 - ► +00. 
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Fig. 7. Shell summed amplitudes of the stream functions shifted by (ih cos 6*, ihsinO*) into 
the most singular direction 6±, for h = 0.0,0.0030,00060,0.0065 (bottom to top). 

To gain additional insight, we now show results in physical space in terms of the real and 
imaginary parts of the vorticity uo = — V 2 F. 8 For this we limit ourselves to n max = 1000. 
As a consequence, k\ runs from —1001 to 0, while k 2 varies from —2002 to 0. We use an FFT 
program with 4096 2 grid points. Fig. 8 shows contours of the real and imaginary parts of 
the vorticity. The symmetries seen are a consequence of dynamically preserved symmetries 
in the initial condition (2). Near the center X\ = ir, x 2 = n there is a highly anisotropic 
large-amplitude sheet-like structure which can be interpreted as the manifestation of a 
smooth singular manifold which gets within roughly four meshes of the real domain. 9 By 
far, the fastest vorticity variation is obtained perpendicularly to the sheet-like structure, 
in the most singular direction 9*. Fig. 9 shows the variation of the vorticity along a cut 
through x± — 7r, x 2 = n in the most singular direction. It is seen that the vorticity becomes 
rather large (around 40 for the real part). The behavior of the real part near the peak, as 
a function of the distance s to the peak is very roughly as l/|s| as seen in Fig. 10. 

It is also of interest to show the variation of the vorticity in the (y±, y2)-plane, that is 
u±( y yi,y 2 ) = uj(tv + iyi,7T + iy 2 )- Symmetry implies that this is a real quantity which, by 
(7), can be written as 

^(y u y 2 )= £ £ (-l)^k 2 F klM e-^ +k ^\ (18) 

k\=— oo k2=— oo 

k — — k-^ ~\~ k 2 . 



8 Since the Fourier transform is real and supported in the quadrant (ki < 0, k 2 < 0) there are 
actually Kramers-Kronig relations between the real and imaginary parts. 

9 This manifold will actually be constructed explicitly in Section 5. 
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Fig. 8. Contours of the vorticity, (a) real part, (b) imaginary part, (c) and (d) enlargements of 
(a) and (b). 

The latter can be viewed as a double power series in the variables W\ = e yi and W2 = e V2 . 
The alternating sign property of the Fk ly k 2 mentioned at the end of Section 3 implies that 
all the coefficients (— l) kl k 2 Fk 1 ,k 2 are nonnegative for ki < —1. It is well known that if 
a power series in one variable with nonnegative coefficients defines an analytic function 
with singularities, then the nearest singularity to the origin is on the real positive axis 
(Vivanti's theorem, see, e.g., Dienes, 1931). Here we expect that there will be a whole 
singular curve in the (y±, ^-plane. The Fourier series can be used to calculate w+(yi, 1/2), 
at least within the disk yf + y% < 5 2 . Fig. 11 shows contours of u+ in and around this 
disk. The contours are almost straight lines perpendicular to the most singular direction, 
because 5 is very small. The variation of uj+ along the the most singular direction as a 
function of the distance y to the origin is shown in Fig. 12. When plotted as a function of 
the distance 5 — y to the nearest singularity (y±+, — 5 (cos 6^, sin#*), it follows roughly 
a power law with an exponent close to —3/4 over almost two decades. We observe that 
the vorticity reaches values around 200. Peak values between 100 and 200 and a scaling 
exponent close to —3/4 are also obtained for the behavior of uo± near other points of the 



11 




Fig. 9. Variation of the vorticity along a cut through (xi, X2) = (vr, it) with the most singular 
direction 0*. (a) real part, (b) imaginary part, (c) and (d) enlargements of (a) and (b). 



singular manifold (to be denned precisely in the next section). We observe that there may 
be an inconsistency between the exponent observed in the real direction which is roughly 
— 1 (Fig. 10) and the exponent —3/4 observed in the imaginary direction (Fig. 12) since the 
vorticity, being analytic, should have the same scaling. The singularities being contained 
in the imaginary plane above X\ = n, x 2 = vr, at least a few meshes away from the real 
domain, the scaling in the imaginary direction is expected to be more reliable. 

Actually, there is strong evidence for the vorticity being infinite at (yi*, 1/2*) coming from 
the behavior of the shell-summed amplitudes. Indeed, we know (i) that all the terms in the 
double sum (18) are nonnegative and (ii) that the shell-summed stream function amplitude 
decreases nearly as k~ 2 over two decades or more (Fig. 7). It follows that the shell-sums 
for the vorticity at (yi+, y 2 *) are almost /c-independent and thus, when we sum over k, we 
get an infinite value. 
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Fig. 10. Same data as in the right half of Fig. 9 (c) in log-log coordinates. To obtain a better 
scaling, we subtract a constant (here 2.0) from the real part of the vorticity. 
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Fig. 11. Contours of the vorticity o>*(yi, 2/2) for x\ = ir and X2 = tt. The dotted circle has radius 
5 = 0.0065, whose smallness makes the contours nearly straight. Contour values from lower left 
to upper right are 30, 40, ... , 90, 100, 200 (the thick solid line is for 100). 

5 The singular manifold 

In one dimension it is well known that the high-wavenumber asymptotics of the Fourier 
transform of an analytic periodic function is governed by the singularities nearest to the 
real domain (see, e.g., Carrier, Krook and Pearson, 1966; Frisch and Morf, 1981). This 
has a non-trivial extension to more than one dimension when the function has singularities 
on a complex manifold S determined by an equation S(z) = 0; let us sketch this in the 
two-dimensional case. Consider a periodic analytic function F(z) given by the Fourier 
series (7). The Fourier coefficients are given by the double integral over the real periodicity 
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Fig. 12. Vorticity uj+ along the most singular direction as a function of the distance 5 — y to the 
singularity nearest to the real domain. 

domain 

Km = tA^ / / dx 1 dx 2 e~ i ^ +k ^F(x). (19) 



Let us set k\ = kcosO and k 2 = ksin.6. We are interested in the behavior of F kl:k2 when 
k — > oo for a given direction defined by the angle 6. For this we change coordinates from 
(x±, x 2 ) to (x\\, x±) chosen, respectively, parallel and perpendicular to k. Note that the 
exponential in (19) involves only x\\. Hence (19) can be rewritten as a one-dimensional 
Fourier transform 



P = 7^j dx\\e- ikx *ge{x{) (20) 
9e( x \\) = ^. J dx ± F(x\\k + x±k±), (21) 

where k = k/k and kj_ is k rotated by +7r/2. It follows from (20) that, for k — > oo, 
F ~ e~ lfcz n where z| is the singularity of ge(<2||) in the complex plane nearest to the real 
domain and with negative imaginary part. 10 Hence |F| ~ e ~ k5 ( e ) where 5(6) = — Im 



Z\\ . 



How do we obtain the singularities of gg(zn), given by the integral (21)? For some (complex) 
z\\ there may be singularities of F(z\\k + z±k±) for real zj_. They can however be avoided by 
shifting the contour of integration away from the real z± axis. If we change z\\ , this will work 
as long as the contour does not get pinched between two (or more) coalescing singularities 



10 For simplicity, we are ignoring algebraic prefactors which depend on the nature of the singu- 
larity. 
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(the pinching, generically, does not take place on the real z± axis). Assuming that the singu- 
lar manifold S can be represented by Sg(z», z±) = in the zu and z±_ variables, a necessary 
condition for this pinching is a double root in z±, i.e. Sg = and dSg/dz± = 0. This system 
of two equations has discrete solutions (zt, z\ ), one of which will control the asymptotics. 
In an earlier version of this paper (available at http://arxiv.org/pdf/nlin.CD/0310044vl) 
we stated that the pinching argument may already be known. Recently W. Pauls pointed 
out to us that it has already been used in a special case by Henri Poincare (see Sections 
94-96 of Poincare, 1899) and has later been generalized by Tsikh (1993). 

All this simplifies considerably for the function F, the solution of the similarity equa- 
tion. Indeed, all the relevant singularities have real part at the point (n, ir) and we can 
restrict everything to the pure imaginary plane (yi, y 2 ) passing through this point. Let 
the restriction 11 of the singular manifold to this plane have the parametric representation 
y±(9), 112(0), where y 1 and y 2 are differentiable functions of the parameter 9, chosen in 
such a way that the angle between the y\ axis and the tangent at the singular manifold is 
9 + 7r/2. We obviously have 

^ = - J-. (22) 
dyi tan 9 



This representation is convenient because the aformentioned condition of having a double 
root is easily shown to express the tangency to the singular manifold of the straight line 
perpendicular to the direction 9 — n at a distance 5(9) from the origin, as illustrated in 
Fig. 13. 12 




Fig. 13. Geometrical determination of the logarithmic decrement 6(0) associated to the direction 
8 in terms of the singular manifold in the (yi, ?/2)-plane. 



The equation for the tangent reads 

5(6) = y 1 cos(9-n) + y 2 sm(9-n). (23) 

11 If, as is likely, the singular manifold is a complex analytic curve, it depends on one complex 
parameter or on two real ones; hence it cannot be entirely in the (y±, 2/2) _ pl ane - 

12 The subtraction of it is because the wavector has negative components and thus lies in the 
third quadrant n < 6 < 3ir/2. 
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Fig. 14. The singular manifold in the (yi, y 2 )-plane constructed as an envelope of numerically 
determined analyticity disks for 305 different choices of parareal planes. Inset: enlargement near 
the origin. 

The system (22)- (23) is easily inverted (provided the singular manifold restricted to the 
(2/1 , ?/ 2 ) -plane is convex) to give 



These equations allow in principle the construction of the singular manifold from the knowl- 
edge of the angular dependence of the logarithmic decrement 5(8), as given in Fig. 6. 

In practice, because 5(8) changes very quickly near its minimum, this is not a very well 
conditioned procedure. An alternative construction by envelope of analyticity disks is now 
proposed. We observe that, when we replace the real plane by a parareal plane shifted by 
an imaginary vector (hi, h 2 ), the location on the singular manifold nearest to this parareal 
plane is generically different from the one nearest to the original real plane. For a general 
2D problem without special symmetry, the singular manifold is a one-dimensional complex 
manifold, and thus may be parametrized with two real variables, e.g. hi and h 2 . For the 
special case considered here, we can work in the (yi, y 2 ) plane and draw around each point 
(hi, h 2 ) an analyticity disk of radius 5 (hi, h 2 ). The latter is determined as the logarithmic 
decrement of the shell-summed average of exp(k 1 hi + k 2 h 2 )F kl;k2 , where the exponential 



V2(8) 



5(8) cosd + 5'(8)sind 
5(8) sm8-5'(8) cosfl, 



(24) 
(25) 



where 5' '(d) 



d5(8)/dd. 
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factor is the consequence in Fourier space of the imaginary translation (ih±, ih 2 ). 13 It is 
easily shown that 



where 5(9) is the logarithmic decrement in the direction 9, defined in Section 4. Since 5(9) 
is rather poorly determined, it is better to measure 5(hi,h,2) directly. We perform this 
calculation for 305 choices of the pair (hi, h 2 ), taken on a regular grid. The shell-summed 
averages are fitted to an exponential with an algebraic prefactor. Only points displaying 
decreasing exponentials, i.e., located below the singular manifold are kept. The result is 
shown in Fig. 14, where the singular manifold emerges as the envelope of all the analyticity 
disks. Let us finally observe that the method of analyticity disks can be generalized to 
problems without any special symmetry; this will be discussed elsewhere. 



6 Concluding remarks 

Let us now summarize what we have learned about singularities for the similarity equa- 
tion (5) governing the short-time asymptotics. First, we remind the reader that for one- 
dimensional problems, complex-space singularities can be at isolated points but need not: 
for nonintegrable ODE's they often form fractal natural boundaries (Chang, Tabor and 
Weiss, 1982). In two dimensions or higher, isolated singularities are ruled out: the singular 
set is a complex manifold (or worse). In the present case, the evidence is that the singular- 
ities are on a smooth and probably analytic one-dimensional complex manifold. This was 
already conjectured in FMB on the basis of numerical results for the full time-dependent 
Euler equation. Now, we have an explicit construction of the singular manifold by the 
method of envelope of analyticity disks (Section 5), which gives strong numerical evidence 
for smoothness. 

The very clean scaling we have observed for the Fourier transform of the solution near the 
singularity (Fig. 7) implies that the vorticity is infinite on the singular manifold Since in 
2D the vorticity is conserved along Lagrangian fluid trajectories (both in the real and the 
complex domain), this result strongly suggests that the singular manifold is mapped to 
complex infinity by the inverse Lagrangian map. 14 



13 These imaginary shifts introduce a bias similar to what is done in the Cramer (1938) derivation 
of the law of large deviations. 

14 Note that this does not imply the absence of complex-space singularities in Lagrangian coordi- 
nates: fluid particles situated initially at finite complex locations can be mapped to infinity later 
on (Pauls and Matsumoto, 2005). 




(26) 
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We have studied here a special flow, whose two-mode initial condition given by (2) has a 
center of symmetry at (jr, ir). Such symmetries are hard to avoid when using a minimal 
number of Fourier modes. At short times, the symmetry constrains the complex-space 
singularity nearest to the real domain to have the real part of its location at (x 1 , x 2 ) = 
(ir, 7r). In FMB we have shown that this ceases to be the case at later times. We do of 
course hope that there is something universal in the nature of the singularities, be it only 
their very existence. We must however stress that in the present case we have a singular 
manifold with a continuous piece and not just isolated points when (x±, x 2 ) = (it, it). This is 
definitely non-generic, since a manifold with complex dimension 1 may be viewed as a two- 
dimensional real manifold in four-dimensional real space. When higher-order harmonics 
are added to the initial condition (2), the situation becomes more complicated and the 
short-time asymptotic equation will in general depend on the ratio 2/2/2/1, that is on the 
direction in which the imaginary coordinates become large. 

Let us finally point out that one particularly challenging problem is to actually prove that 
there are singularities in the complex domain: at the moment we are not aware of a single 
instance of a solution to the incompressible Euler equation with smooth initial data for 
which the existence of a singularity (real or complex) is demonstrated. In the present case, 
the proof may be facilitated by the observation made in Section 4 that, in the (2/1, 2/2)- 
plane the vorticity is real and can be represented as a double power series with (apparently) 
nonnegative coefficients. A lower bound on these coefficient - which are obtained from the 
solutions of the recurrence relations (12) - could lead to such a proof. 
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